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The propagator of confined quarks is calculated for timelike momenta by transforming Minkowski 
Greens functions to the Temporal Euclidean space. Based on the framework of the Schwinger- 
Dyson equations the QCD quark propagator is obtained in two approximations which differ by 
assuming behaviour of gluon propagator. In both studied cases we get universal result for the light 
quarks: The quark mass function becomes complex bellow expected perturbative threshold, the 
obtained absolute value of the infrared mass is M ~ Aqcd with the infrared phase ~ 5 . Permanent 
confinement of quarks is maintained by generation of the complex mass function which prevents 
■ a real pole in the propagator. We will show that timelike dynamical Chiral Symmetry Breaking 

(CSB) solution is approximately, but non-trivially determined by the solution of gap equation in 
f^) the standard Euclidean space. 
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I. INTRODUCTION 



Quark confinement in Quantum Chromodynamics (QCD) is theoretically unsolved phenomenon of longstanding 
interest 0, 0, 0, B S S H 0|- QCD as the strong coupling theory is not easily tractable, thus the confinement of 
colored object has been questioned in various indirect ways. Considering static infinitely heavy color sources the 
linearly rising potential was predicted and observed by simulation on the Euclidean lattice [l(| • 
i' The alternative to the lattice theory is the continuous framework of the Schwinger-Dyson equations (SDEs) which 
in principle could provide a unique powerful tool for the nonperturbative QCD study (for a reviews see [3 El El, 
4^ i ' EH)- While it is naturally expected that the SDEs Euclidean formulation should approximately corresponds to the 
continuum limit of lattice data whenever they are available, however the real advantage of SDEs approach is the 
possibility to explore directly all Minkowski space. To that point, based on the spectral representation, the problem 
has been formulated in Minkowski space and the results have been actually obtained for a weak or a medium coupling 
quantum field theory (for a review see [IB]). Actually for several cases of toy models the true equivalence between 
Minkowski Greens functions evaluated at spacelike momenta and the calculated in the Euclidean space form the 
beginning has been explicitly checked by the numerical solutions. For gluodynamics (QCD without quarks) the 
problem was formulated by using covariant Gauge/Pinch technique method time ago [la ], however never solved in 
Minkowski space in practise. Recently, we have no idea how much the analytical assumptions are justified for gluons, 
however we have no doubts that the usual analyticity is too strong assumption for the quark propagator [l7j . To get 
chiral symmetry breaking with an analytical quark propagator, i.e. which has a standard branch cut on the positive 
• • ■ p 2 semi-axis, is very likely impossible - the quark propagator doesn't satisfy Lehmann representation. 

Direct solution of SDEs for strong coupling theory system in Minkowski space is still very complicated task and 
as we discuss above it requires an additional assumptions. To get rid of part of the mentioned weaknesses (e.g. 
}_i ' assumption of the spectral representation) the SDEs has been recently formulated in Temporal Euclidean (ET) space 
[l8l [l3]. providing us with the solution for Greens function for the timelike, but Euclidean momenta. In this case we 
have the Euclidean metric p 2 ET = Pi 

P= (pi,P2,P3,Pi) [ET] = {po,-Wi,-ip2,-iP3) [M] , (1-1) 

where Eq. (jl.ll ) explicitly shows the prescription between the components of Temporal Euclidean -ET- and Minkowski 
-M- fourvector. Recall here, the prescription for the measure d A pM — id ET with real ±oo boundaries in a momentum 
loop integrals. Likewise for the standard -time component- Wick rotation [2(| the equivalence between ET space and 
Minkowski subspace is the assumption. Here, it is noteworthy that it was recently proved in fl9j ] that the ladder 
approximation of Minkowski QED2+1 is exactly equivalent to the formulation of the problem in ET space. 

If the full fcrmion propagator has no mass singularity in the timelike region, it can never be on-shell and thus never 
observed as a free particle [2ll . |23| . |24| . [29( . In our solution the imaginary part of the mass function is automatically 
generated for a coupling strong enough. In the paper [l8[ such result was firstly obtained for explicitly massive quark 
with bare current mass m ~ Aqcd- At that time the authors of [l|| were not able to obtain the results for an 
arbitrary quark mass because of numerical obstacles accompanying their specific model. Therefore the main purpose 
of our paper is to present a models and techniques that exhibit good numerical stability for all quarks flavors, e.g. for 
the light u, d quarks and hypothetical massless case as well. We will show that the light quark mass is the complex 



2 



functions with relatively tiny real parts in the infrared. The infrared complex phase is approximately given by the 
ratio of the current and the infrared (constituent) quark mass. Quite interestingly, described complcxification is 
universal and it doesn't depend on the details of the interaction kernel. As we will show explicitly the qualitative 
feature of observed complexification is independent on the (un)presence of (pole) type singularity in the kernel. 

As an another interesting model possibility we decreased the infrared running coupling and we study confinement 
and CBS phenomena near the phase transition, where generated mass is particularly small when compare to the scale 
of confinement, M << A. The observed results are relevant for the Technicolor models, eg. to walking Technicolors 
wherein large number of fermions makes the strong coupling softer in the infrared (originally these models have been 
developed to avoid flavor changing neutral current, however they have theirs own interest). 

After the introduction of the models details in the next Section the numerical results are presented in the Section 
III. In this section we discuss the relation with the standard Euclidean solution. In the Section IV we summarize and 
conclude. 



II. THE MODELS- DRESSED LADDER APPROXIMATION OF QUARK SDE 

The quark gap equation we solved is the rainbow ladder approximation of the full quark SDE. In this Section we 
describe the details of the models employed here. In the ladder approximation the quark SDE reads 

S(p) = *C A g 2 J -0^ 7a G^(p- q )S(q) l0 , (2.1) 

where the full gluon-quark-antiquark vertex of the exact SDE has been replaced by the standard four dimensional Dirac 
matrix, which obey anticomutation relation {7^, 7^} = 2g M „ (we use the Minkowski metric g M „ = diag(l, —1,-1,-1)) 
and where C A — T a T a = 4/3 for SU(3) group and G a/3 is the gluon propagator. The quark propagator S can be 
parametrized by two scalar functions conventionally like 

Sip) = S v (p)^ + S.( P )= M(p) 1 _ B(py (2.2) 

The renormgroup invariant mass function is defined as M = B/A, the physical mass M p could be identified by the 
pole position of S ,i.e. by the solution M 2 (M 2 ) = M 2 in unconfining theory. 

We assume these functions are complex, therefore it is convenient to parametrized them in the following way 
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where R A , Rb (T a ,Tb) are the real (imaginary) parts of the functions A, B and the denominator D reads 

D = {[R\ - T\]k 2 - [R% - T|]) 2 + A(T A R A - T B B) 2 . (2.5) 

The last missing ingredient which completes our considered gap equation is the gluon propagator. With the 
exception of the large timelike momenta p 2 > Aq C£) , where the effective coupling is small and the result is available 
by the analytical continuation of the Euclidean perturbation QCD, gluon propagator is basically unknown function 
for timelike momenta. Most of the last decades nonperturbative studies were devoted to the quark SDE in Landau 
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gauge, thus for a possible comparison we will work in this gauge as well. The Landau gauge Gluon propagator is 
completely transverse 

G> = ~ g ^ + ( *~°^ G(fc 2 ) , (2.6) 
and is fully determined by the gluon form factor G which we will model as it is described bellow. 



A. Model I. 



The model I. is based on a simple generalization of the perturbative one loop result for the gluon form factor. 
For this purpose we consider the kernel (in fact the effective product of gluon propagator and 7^ part of the vertex 
function) such that it has a standard single pole and the the function G is taken as 



47T 



G(q 2 ) 



4tt/ /3 







\ln 





(2.7) 



where in (|2.1ip represents the beta function coefficient, for which we take Air/ (3 = 1 (recall, 4ir/f3 ~ 1.396 for three 
active quarks in perturbative QCD). The prefactor is adjusted in a way that Eq. (|2.7p behaves as QCD running 
coupling at ultraviolet, i.e. 



» Aqcd) 



4tt//3 



ln(±q*/A* QOD ) 



(2.8) 



where ± stands for timclikc or spacclikc q 2 respectively. In the infrared effective running coupling defined as 
g 2 G(Q)/ (4tt) = 2, thus the function G is large enough to generate QCD typical CSB and and it provides confinement 
of quarks for any flavor as well. 



B. Model II. 

The second model we use represents a slight modification of the one already considered in [l8[ . In this case, the 
gluon propagator is expressed through the following spectral representation: 

^ = r^, (2 . 9 ) 

47r q l J q A - v 

The function p g we consider here is a regular smooth function, ensuring thus there is no pole in the gluon propagator. 
The main difference comparing to (l8j is that we drop out the usual Feynman ie prescription so the principal value 
integration is understand for the timelike momenta in (|2.9[) • As far as we are not considering the feedback of complex 
quark loop the gluon propagator is clearly real for all Minkowski q 2 . 

In principle the weight function p g can be obtain by solving the gluon gap equation as suggested in [l(| and further 
considered recently in (2f| [26|, [27| (up to the changes followed from the ie absence). 

Following the arguments presented in [2?| . it seems that gauge invariant gluon propagator is not "strong" enough 
to trigger CSB. Very likely, the next leading vertex correction must enforce the kernel to get the expected picture 
of CSB and confinement. To that point we simplify and make a phenomenological choice of the function p g which 
provide the desired solution. 

The function we actually use in our numerical study reads 

a{x) p a (x) 

^^MojxTolA 2 (2 - 10) 



where the function a(x) is calculated through 

Pa{x) 



n 2 -ln 2 (x/A 2 QCD ) 



a(x) = P. I dv 9 -^, (2.11) 







4 



where symbol P. stands for Cauchy principal value integration. For the reader familiar with our previous paper on 
Temporal Euclidean QCD, now gluon form factor is twice weaker there is factor 2 omitted in (|2.11[) . further the 
infrared cutoff 0.1A 2 is introduced in (|2.10| in order to have finite and smoothed gluon propagator G/q 2 in the 
infrared. Evaluating the integral one can see that for a very large momenta the gluon propagator is softened by the 
power of log, i.e. 



4n {lQ 1 A Q cd) ~ 2ln 2 (±q 2 /A 2 QCD ) 



G(\q 2 \ » A 2 QCD ) ~ ,^J^ 2 , (2.12) 



(this fact has been overlooked in pil). This affects the UV tail of quark mass function, but is quite unimportant for 
the low the energy behaviour of the quark mass function M(0). 



III. SOLUTION OF THE GAP EQUATION 
A. QCD light quarks 

The quark gap equation have been solved numerically by the standard method of iterations. The Gaussian numerical 
integrator with 600 mesh points was used for all the numerical integrations. The angular integration has been 
performed analytically in the model II, thus we only numerically integrate over the auxiliary variable v. To avoid 
numerical noise, large density of the mesh points for q 2 ~ A, v ~ A was used in the Model II. Furthermore, in order to 
achieve numerically stability during the run of iteration process, we gradually decreased current quark mass to reach 
the massless limit. The model I is completely stable, whilst there is a small numerical noise for a low q 2 in the model 
II. 

Actually, what we have solved was the system of four real equations for the functions Ra, Ta, Rb,^b- These have 
been obtained by standard trace projections (of Dirac matrices), what leads to the coupled system of two equations 
for the functions A, B in the first step.. After performing the 3d-space Wick rotation, we get the expressions in the 
Temporal Euclidean space. Further using Eqs. (|2.4p we arrive to the equations for the imaginary and real parts - 
Cartesian complex coordinates of A, B - the functions Ra, ^a, Rb,^b ■ Procedure is very straightforward and we do 
not list the result here, however the reader can find the equation for B in Section [ill CI where the properties of the 
solution are discussed. 

As in the case of perturbation theory, the equations contain UV divergences which require renormalization. The 
chiral limit is considered for both the models. The only renormalization function A is renormalized, the function B 
is finite in this case. In the case of light quarks we assume that the ultraviolet timelike behaviour is similar to the 
spacelike one, hence we set a small real renormalized quark mass M(/i 2 ) ~ Aqcd/W0 for a very large timelike /i 2 . In 
this paper all dimensionfull quantities are scaled by Aqcd, having a typical value for Aqcd = 250 — ZhQGeV we set 
a few MeV light quarks at 200 — 300GeV (a given mesh point is chosen for convenience, thus the corresponding value 
of the renormalization point is /i = 920.3Aqct> f° r presented solution here). In this case we subtract the real parts 
of functions A, B setting the renormalized values ReA — 1, ReB — Aqcd/100 at the renormalization scale fi. As in 
[l8|, the imaginary part is expected to be finite and therefore not subtracted. This procedure clearly maintains the 
hcrmicity of the classical Lagrangian. 

Having the numerical solutions we calculated the magnitude and the mass function phase defined as M = B/A, M = 
\M\e l ^ M . The appropriate results are presented in Fig. 1-3. Renormalization wave function A for our models is shown 
in Fig. 4. It is approximately real everywhere, the imaginary part portion obtained ~ 10 -8 is smaller then estimated 
numerical error. There is no remarkable difference between the solution of A for u, d quarks and for the exact chiral 
limit. In the case of model II these two lines are not distinguishable. The absolute values of the mass functions are 
shown in Fig.l. In the case of the model II, the function \M\ shows up the maximum at ~ 2Aqcd where it also cuts 
the linear function of p. This details are better seen in the Fig. 2 where we show the infrared details with linear axis 
scaling. 

As a bonus of our CSB solution we get the pion Bethe-Salpeter wave function %(P = 0,p) ~ B(p) [ll[, which is 
the Goldstone boson manifestation of broken chiral symmetry here. 

Not similarly to previously studied confining theory QED2+1 [l9l |. where a quite tiny imaginary part preserves 
regularity of the time axis propagator, here a huge dynamical generation of imaginary part mass function M is 
observed. For pure CSB (to = 0) case the observed generated mass is purely imaginary and the same is approximately 
valid for the infrared mass when small explicit breaking is considered. In both cases, the imaginary part vanish at 
high p. The phase <f> of the running quark mass function for the models I and II respectively is shown in Fig [TJ The 
chiral limits are always <f> = 7r/2 for both models and hence not displayed. Increasing a current quark mass, the phase 
is decreasing (rather say its absolute value, there are possibly more than one solution for a quark light enough [28|), 
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FIG. 1: Magnitude \M\ of the running quark mass function M = |M|e 10 for modeled QCD 1,11, its chiral limit CI.CII. The 
"Walking Technicolor" T solution is added for for the comparison, scale is Aqcd = 1 (and Arech. = 1 is set up in this case as 
well) 

however from a certain value of the quark current mass it does not disappear at its stay constant from some p. A 
more complete study on heavy quark confinement will be published elsewhere. Further discussion of the solution, 
after the following Technicolor digression, will follow in the section IIII CI 



B. Criticality in large Nf QCD, walking Technicolor and views in ET space 

QCD is an example of non-Abelian gauge theory with small -two or three- approximately massless fermions. If the 
number of massless fermions is larger, but the asymptotic freedom is still preserved at some high scale, the theor y is 
conformal in the infrared, which have been already suggested by the analysis based on two loop beta function (3fj|. l3l| 
and analyzed for variety of gauge theories (see [32| and references therein). Non-Abelian quantum filed theory with 
possibly large but yet supercritical number of (classically) massless fermion appears as an alternative candidate for 
electroweak symmetry breaking [33l. |34L |35| . Typically in these "walking Technicolor" models, the number of fermions 
set up the scale At which characterize running of the effective coupling, which is much larger then fermion mass 
generation M(0) < Ay. In other words, the conformal infrared nontrivial coupling is adjusted a few percentage above 
its critical level , say denned at zero momenta a c (0) = g 2 /(At:) , where g is a gauge coupling. In these theories the 
coupling is "slowly walking" with momentum in the infrared p < A where it ensures CSB. The CSB and Technihadrons 
spectrum in these models have been studied in (36l . l37l . l38l . |39| . 
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FIG. 2: I infrared behaviour of the functions M as they are in the Fig. 1, but in with linear axis. The solution for large Nf is 
omitted here. 

Identification of the critical number of the flavors for given theory requires nonperturbative knowledge of infrared 
behaviour of the running gauge coupling. For SU (3) gauge group a recent estimate based on the lattice simulations 
gives 8 < Nf < 12 l4l| To model the quark gap equation for large Nf we adjust the infrared gluon form factor 
to be close to minimal strength necessary to trigger CSB. It is achieved by taking the constant = 1/3 of the model 
/, which model is particularly suited for this purpose. We have found the critical couplings are the same in the 
Standard Euclidean formulation and the one obtained in the Temporal Euclidean space. The obtained CSB solution 
in ET space solution is confining one, the mass function B is purely imaginary in his case. The renormalization wave 
function A remains real receiving expected smaller corrections than in the case of QCD. 

Decreasing infrared coupling furthermore we obtain only trivial solution for the function B. 

C. Further observation 

Let as consider the gap equation in more details. As the consequence of our special model I we get with astonishing 
accuracy the following relation 

A( P 2 > 0) = A ET ( P 2 ) = A E ( P 2 ) = A( P 2 < 0) 

B(p 2 > 0) = B ET (p 2 ) = iB E {p 2 ) = iB(p 2 < 0). (3.1) 
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FIG. 3: Phase (f> of the running quark mass function M = \M\e 1 ^ for the models I and II respectively, axis momentum is in the 
units of Aqcd- 



between spacelike and timelike solution. The first is obtained in standard Euclidean formulation, while the second 
in ET space for pure CSB, m = 0. A small current mass leads to a small deviation only. 

The exact symmetry (|3.ip is directly visible by the inspection of the SDE for B which in ET space reads 



4 t~ 
B(x) = - \ 

T Jo 



dy 



B{y) 



A 2 (y)y-B 2 (y) 



dzV{x, y, z) 



V(x,y,z) 



yVl - z 2 



q 2 log(e + g 4 /A 4 ) ' 



(3.2) 



where here q 2 = x + y — 2^/xyz. 

For the standard (spacelike) Euclidean formulation the gap equation reads 



4 f°° 
B(x) = - dy 
it Jo 



B{y) 



A 2 (y)y + B 2 (y) 



dzV(x, y, z) 



(3.3) 



The symmetry 13. II is manifest. In other words: for the Minkowski kernel which is even with respect of argument of 
the gluon propagator G(q 2 ) = G{— q 2 ), the propagator functions Ss obtained in E and ET spaces differ by the phase 
factor 



5 s (V)=e ±W2 5 s (p 2 ), 



(3.4) 
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FIG. 4: Renormalization wave functions for the model I and II. The same is shown for the chiral limit (which are indistinguishable 
in one case). The "Technicolor" T solution is added for for the comparison, the momentum axis is scaled in Aqcd and At 
respectively. 

(we use Minkowski space convention in this expression) while the the functions SV's are identical in both spaces. 

The gap equations above are displayed for small (3 = 1/2 and we confirm the solution numerically for any beta, i.e. 
the symmetry is kept for for increasing /?, when we are gradually leaving the critical point and reaching the value 
typical for QCD. No other solution was observed numerically for purely dynamical CSB. We observed that the non 
zero real imaginary parts of B and A are generated because of nonzero current quark masses, therefore the observed 
identity (|3.1I are only approximate there in the infrared. In reality, the full gluon propagator and the quark-gluon 
vertex are not suppose to be a simple real and even functions of momenta, but we expect substantial changes due to 
this. We argue here, solutions of quark gap equations for small current masses already presented in the vast amount 
of the literature are not only Euclidean solutions usually identified with spacelike Minkowski solution, but up to the 
phase, they represent a rough but reliable approximation of the timclike Minkowski solution as well. 

Although we assume here, the observed structure is crucial for selfconsistent determination of the Greens function 
in ET space, the resulting phases become irrelevant whenever observables are composed, since the scalar product 
made from a given amplitude and its conjugated can only represents a measurable quantity. 
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IV. SUMMARY AND CONCLUSIONS 

The first analysis of the light quark gap equations in the Temporal Euclidean space was presented in the paper. We 
assume that ET solutions represent Minkowski solutions at timelike axis as good as the standard Euclidean solutions 
the Minkowski spacelike one. Behaviour of the solutions were discussed, the main feature of the all solutions is that 
the quark mass function is predominantly imaginary, providing confining solution for the quark propagator; such 
quark propagator has no pole nor branch point at real p 2 axis. 

In perturbative QCD and in weak coupling theories generally, the well known calculational trick- the Wick rotation- 
is under good control. It is a text book knowledge that timelike solution could be obtainable by an analytical 
continuation of the result conveniently defined and calculated in spacetime Euclidean space. Here, we have assumed 
no poles and no cuts at a real axis and got complex solution justifying our assumptions. The observed CSB solution 
in ET space cannot be obtained by analytical continuation of the Euclidean one, it basically differs by the phase. The 
meaning of that requires more deep understanding. 

Recently we do not know how to judge and evaluate the quality of an assumptions we made when switching 
between Minkowski and Euclideans worlds. The all complex space of four momenta is not under easy control. On the 
other side, the Euclidean metric simplifies things, how well it approximates our real Minkowski world can be judge 
a posterior: The hadrons properties calculated from the QCD Green'S functions are recent (future) prospectors of 
quality of Euclidean (Temporal Euclidean) calculations. 

Some possible doubt can follow from unknown prescription of the kernel of considered SDE here. The gluon SDE 
and related vertices require future selfconsistent analyzes in ET space. It will be certainly more evolving but similar 
task as the one performed here for the quark gap equation alone. On the other side, we do not think that the 
improvement towards the exact knowledge of the SDE kernel can drastically change qualitative result. The observed 
complexification phenomena is very universal as it has happened to two quite different models. However, going beyond 
the rainbow ladder is task for future Temporal Euclidean space study. 

At last but not at least the dynamical CSB has been studied near the critical coupling. The observed numerical 
solutions obtained do not suggest separation of the dynamical CSB from confinement. These phenomena go hand by 
hand in hypothetical Walking Technicolors models. 
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